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o I Optical soliton pulses offer many applications within optical communication 

systems, but by definition a soliton is only subjected to second-order anoma- 

^ lous group-velocity-dispersion; an understanding of higher-order dispersion is 

^ necessary for practical implementation of soliton pulses. A numerical model 

GN of a waveguide was developed using the Nonlinear Schrodinger Equation, 

^ with parameters set to ensure the input pulse energy would be equal to 

^ the fundamental soliton energy. Higher-order group-velocity-dispersion was 

^ gradually increased, for various temporal widths and waveguide dispersions. 

• • 

>► A minimum pulse duration of 100-fs was determined to be necessary for fun- 

^ damental soliton pulse propagation in practical photonic crystal waveguides. 
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1. Introduction 

The purpose of this effort is to, both theoretically and with modeling and simulations, study 
the effects of higher-order dispersion on fundamental soliton pulse propagation of ultrafast 
pulses propagating within an optical waveguide. Optical solitons are undistorted standing 
waves that can propagate over long distances undistorted [1,2], which result as part of the 
interplay between the nonlinear Kerr effect and anomalous dispersion [1-7]. These pulses are 
of interest to the scientific community for their self-sustaining properties, which are useful 
both for low-powered optical data transfer as well as high-energy pulsed laser systems. These 
pulses have been generated in optical fibers [1,6], where the stable nature of the pulse reduces 
loss and cross-talk. This effort is to study the feasibility of self-sustaining soliton pulses in 
the presence of higher-order group- velocity dispersion (GVD) [1], with emphasis placed on 
practical studies of photonic crystals waveguides (PhCWG) [3,8-11], which are subjected to 
much higher-order dispersion as a result of the geometric design of the different waveguide 
dielectrics. 

The optimal way of making such a determination, other than through experimentation, is 
to use numerical methods. These methods include Finite Difference Time Domain (FDTD) 
and the split-step Nonlinear Schrodinger Equation (NLSE) method [1,12-15]. The NLSE is 
one of the most effective methods for numerically determining soliton pulse evolution in a 
nonlinear dispersive waveguide, as it is easy to implement and computationally efficient. Just 
like FDTD, the NLSE method is derived from Maxwell's equation [2], under the assumption 
that the Slowly Varying Envelope Approximation (SVEA) is applicable [2, 16], where it is 
assumed that the envelope of the pulse varies slowly in time and space compared to an 
optical period. 

Soliton pulses within fiber optics is a well-established subject, and it has been used pre- 
viously for data communication over many kilometers of fiber networks [1,6]. The primary 
disadvantage of fibers is simply the long-length scales required for the optical nonlinearities 
to occur. There would be many advantages to having soliton pulse compression occur at 
the millimeter length scales, which would allow for better optical data transfer and process- 
ing within a typical silicon computer chip [3-5, 12, 15, 17-22]. For this reason, there has been 
great interest in the development of photonic crystals, which possess extraordinary nonlinear 
characteristics at much smaller length scales. 

One of the first examples of soliton pulse compression at the millimeter length scale has 
been observed in gallium indium phosphate (GalnP) PhCWG [3]. Ultrafast pulses measured 
at 2 picoseconds (ps) were propagated through these PhCWG, which are engineered to 
have extremely high anomalous dispersion [7], and the output pulse was measured with 
autocorrelation at 500 femtoseconds (fs). These PhCWG, by their definition, were measured 
to have orders of magnitude greater second and 3^^-order GVD coefficients when compared 
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to silica fibers [3,9]. 

One issue inherently manifested in photonic crystals is the robustness of NLSE and the 
soliton wave; due to the highly dispersive nature of photonic crystals the SVEA tends to 
break down quickly. The robustness of the NLSE for few-optical cycle pulses has been studied 
both experimentally and numerically in optical fibers [16,23], and analytically in photonic 
crystals [10]. It is clear that proper engineering of the higher-order dispersion in the PhCWG 
is necessary for sustaining the shape of the pulse as it propagates through the waveguide. 
To the authors' knowledge, however, there has been no formal study to observe when an 
ultrafast pulse is too short for the soliton wave to sustain itself in a highly dispersive waveg- 
uide. For this reason, this study was conducted, where a fundamental soliton was propagated 
through a highly dispersive waveguide, and the sustainability was studied for various tem- 
poral durations and higher-order GVD. 

2. Basic Theory 

In order to accurately simulate the optical soliton propagation, two crucial nonlinear char- 
acteristics are necessary; the GVD and the nonlinear Kerr effect [1,2]. The GVD coefficients 
are determined by knowing the derivatives of the refractive index over the angular frequency. 
These terms are found using the following equations [1]: 

^{u) = n{u)-, (1) 

Co 
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where Cq is the speed of light, uo is the angular frequency, vq is the group velocity, and uq is 
the group index. These terms can be determined either experimentally [8,24], approximated 
with the Sellmeier equation [1,2], or engineered to different values as is the case with photonic 
crystals [9]. 

The other component of generating soliton pulse compression is the nonlinear Kerr effect. 
This results in a change of refractive index proportional to the square of the electric field, 
or proportional to the electric field intensity. The change in refractive index is determined 
by [1,2,17]: 

An = n2/i = h, (5) 

eon 
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where n2 is the nonhnear Kerr coefficient (m^/W), Ii is the intensity of the pulse (W/m^), 
n is the refractive index, eo is the electric permittivity in a vacuum (8.854 • 10~^^ F/m), 
r]o is the vacuum impedance (r^ 120-7r), and x*^^) is the 3'""^-order nonlinear susceptibility 
coefficient [2]. These two nonlinear parameters can combine together to form the NLSE: 



7 = ^-(^r, (7) 

where (s^/m) is the 2'''^-order GVD coefficient, /3s (s^/m) is the S'^'^-order GVD coefficient, 
Ae (m^) is the effective area of the waveguide, A is the wavelength (m), 7 (m-W)~^ is the 
nonhnear parameter, and a (m~^) is the hnear loss coefficient (disregarded throughout this 
theoretical study). 

When the pulse is focused by its own intensity, it is known as Self-phase Modulation 
(SPM), and when there is strong SPM interplaying with strong anomalous GVD, then soliton 
pulse propagation may occur [3] . The SPM will cause the leading edge (longer wavelengths) 
frequencies to be lowered, and the trailing edge (shorter wavelengths) frequencies to be raised. 
At the same time, the anomalous GVD will cause the lowered leading edge frequencies to 
slow down, and the raised trailing edge frequencies to speed up. This results in the pulse 
narrowing, or compressing into what is known as the soliton pulse [1-3, 19]. 

3. Numerical Modelling 

The NLSE method is also known as the split-step method, as each incremental distance step 
involves two fundamental steps. The first step is to account for the GVD. The code will, 
using the Fast Fourier Transforms (FFT) technique, convert the pulse information into the 
spectral domain, apply the GVD, and then convert the pulse back from the spectral domain 
to the temporal domain. The equation for the GVD is as follows [1,2, 14,25,26]: 

A{z,uj) = A{0,uj)'exp[i^/32^'^z + i^/3suj^z + ...], (8) 

/oc 
A{z,T) ■exp[-2mujT] dT, (9) 
-OO 

T ^ t--, (10) 



where A{z,T) and A{z,<jo) are the pulse envelope functions in the temporal and spectral 
domains, T is the normalized time, Vg is the group velocity, z is the propagation distance, 
and oj is the normalized angular frequency. The second step is to convert this dispersed 
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pulse, convert it to the temporal domain, and apply the effects of SPM. The equation for 
the nonlinear phase shift is [1,2, 14]: 



A{z,T) = A{0,Tyexp[i{—ynr\A{0,T)f]. 

^0 



(11) 



By following through these two steps with each distance increment, one can predict the non- 
linear effects on a soliton pulse envelope during and after optical propagation in a waveguide. 

The first step in this effort was to develop a numerical model to solve the NLSE. What 
was unique about the model used in this study was that it adjusted the input pulse energy 
to consistently propagate a fundamental soliton, depending on the related parameters such 
as the 2^^-order GVD, the Kerr coefficient, effective area, wavelength, group index, effective 
index, and input temporal duration. The input power was calculated by the equation [1,3]: 



■T 



(12) 



where il^ is a constant to account for the numerical error, and r is the temporal duration in 
seconds. As expected, it was observed that for a given soliton number, N'^ = Pulse Energy / 
Fundamental Soliton Energy [3], the pulse shape would look identical after a given number 
of nonlinear lengths. 




Pulse Duration 





Pulse Duration (ps) 



Propagation Distance (mm) 



Fig. 1. Ideal Soliton Propagation, for the (a) fundamental soliton N = 1, (b) 
higher order soliton N = 2, and (c) higher order soliton N = 3. 

As expected, the self-focusing from the Kerr nonlinearity focused the pulse, and canceled 
out the anomalous dispersion of the waveguide. In the case of the fundamental soliton, the 
pulse experienced little change over distance, whereas with second and third order solitons 
there was compression followed by splitting, with the pulse repeating this phenomena while 
propagating through the waveguide. The results were verified to be the result of soliton 
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propagation, as the pulse broadened dramatically when the Kerr coefficient was turned off 
as the pulse propagated in the waveguide. 

An interesting property of the fundamental soliton that was observed numerically was the 
fact that, for a specified soliton number and nonlinear length, changing any of the important 
nonlinear parameters, such as the 2^^-order dispersion, the Kerr coefficient, effective area, 
wavelength, group index, effective index, or temporal duration, the pulse would experience 
no change in shape; only the scales would change. This is assuming there are no losses, 
whether linear or nonlinear such as two-photon and free-carrier absorption [3, 15], as well as 
no higher order dispersions. While the effects of 2^^-order GVD will inherently be greater 
with shorter pulses (which mathematical must have larger spectral widths [27]), by adjusting 
the energy accordingly to maintain the fundamental soliton, the pulse may continue to be 
self-sustaining. The model demonstrated, based on the theory of NLSE, that a soliton can 
propagate through a waveguide, regardless of the waveguide length, dispersion, or temporal 
duration. 

4. Simulations of Higher-Order Dispersion 

The first step in this effort is to focus on higher-order nonlinearities as they disturb funda- 
mental soliton pulse propagation in an optical waveguide. The model was set-up in order 
to determine the changes in Time Bandwidth Product (TBP) [28] of fundamental solitons 
propagating with a range of 2^^-order GVD coefficients, spanning from low-dispersion fiber 
waveguides to highly-dispersive PhCWG. The study encompassed 3^^, 4^^, 5^^, and 6*^-order 
GVD, in order to determine the maximum higher-order GVD term that the fundamental 
soliton can sustain itself over a given nonlinear length scale. 

The simulation was built to analyze waveguides with 2^^-order GVD coefficients that 
ranged from 10~^ to 10 ps^/mm. As this simulation is purely theoretical, the model only 
varied the orders of magnitude. The simulation ran in 2 fs increments from 2 fs to 1 ps. 
Each higher-order GVD coefficient was increased independently and incrementally from no 
higher-order dispersion to 0.1 ps^/mm, 0.05 ps^/mm, -0.3 ps^/mm, and -1 ps^/mm, as the 
maximum 3^^, 4*^, 5*^, and 6*^-order GVD coefficient. Equivalent waveguide lengths of 0.01 
to 5 nonlinear length ratios were used in the study. 

After the simulations, the data was analyzed to determine the maximum higher-order 
dispersion that could be tolerated for a fundamental soliton of a given temporal duration. 
At each specified pulse duration, nonlinear length, and 2^^-order GVD, the fundamental 
soliton was considered to be disturbed if the TBP changed by more than 10%. It was 
observed that the waveguides with highly-dispersive 2^^-order GVD values were more 
tolerant to a given higher-order GVD. This is expected, as a fundamental soliton within 



6 



a highly-dispersive waveguide will have a greater input energy and thus greater SPM to 
balance out the increased group-velocity dispersion. In addition, longer waveguide lengths 
were observed to be more sensitive to higher-order GVD; this is expected as the dispersion 
is proportional to the propagation length. 

After the simulations were conducted, an effort was made to determine if a mathematical 
function could be found of the maximum higher-order dispersion as a function of the 
temporal width, nonlinear length, and 2^^-order GVD. This is a straightforward linear 
algebra problem utilizing least squares [29] : 

Ax = b, (14) 
A^Ax = An, (15) 

{A^Ay\A^A)x = x = (i^i)"'i^6, (16) 

where A represents the values of the individual functions to curve fit with the independent 
variables (temporal width, 2^^-order GVD, and nonlinear length ratio), is the transpose of 
A, X represents the constant coefficients of these functions, and h represents the dependent 
data of the study (the maximum higher-order GVD). The function found to most-closely 
represent the data is: 

/(r, /32, L^l) = Cr/i(r, /32, L^l) + Crf2{T, /?2, L^l) + C3, (17) 

/(r, ^2, Lnl) = Cvr-prexpi-LNL) + C2-T-^2-exp(-2-LjvL) + C3. (18) 





3'''^-order 


4*''-order 


5*'^-order 


gtii_order 




-2.712-10-1 


-2.186-10-1" 


1.277-10-25 


4.898-10-37 


C2 


2.198-10-1 


1.776-10-1" 


-1.006-10-2^ 


-3.939-10-37 




1.377-10-3^ 


2.112-10-"^ 


-1.138-10-5^ 


-3.236-10-™ 



Table 1. Coefficients for Equation 18 



By solving the least-squares equation, one can find the coefficients of the chosen equation, 
which are listed in Table 1. The average error between equation 18 and the simulated 
maximum higher-order GVD tolerances was found to be less than 3%, thus validating the 
equation. By following equation 18, one can reasonably predict the maximum higher-order 
dispersion possible for a fundamental soliton propagating in a waveguide with a given 
2^^-order GVD, temporal width, and nonlinear length ration. 
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(a) 

Minimum 3rd Order GVD for fundamental soliton propagation 



(b) 

iVIinimum 4th Order GVD for fundamental soliton propagation 




Fig. 2. Maximum (a) 3'^, (b) 4^^, (c) 5^^, and (d) 6^^-order GVD coefficients 
for a 10 fs pulse, for a given nonlinear length and 2^^-order GVD. 



5. Photonic Crystal Waveguide Analysis 

In practice, there will always be some higher-order dispersion, and in a PhCWG this can be 
significant enough to measure experimentally. In the analysis of the GalnP chip, the third 
order dispersion (TOD) was measured up to 0.2 ps^/mm [3]. The next step in this effort 
was to run numerous NLSE simulations of the fundamental soliton, for various temporal 
durations and TOD, and to determine whether or not the self-sustaining characteristics can 
exist in the presence of TOD. 

The model was first set at ten nonlinear lengths. As expected, when there was zero TOD, 
there was no change in output TBP or output phase standard deviation, regardless of the 
input pulse duration; the TBP of the output was always equal to that of the input. When 
the dispersion was increased to 0.1 and 0.2 ps^/mm, there was a noticeable change in the 
temporal duration for input temporal durations less than 100 fs. The TBP would vary from 
half the input TBP, which is evidence of higher-order soliton compression as a result of the 
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increased dispersion; and up to over twice the TBP, which is evidence of even greater pulse 
broadening and distortion of the pulse shape. 



tl. 

Pulse Duration (s) 

Fig. 3. The output TBP over the input TBP for an optical pulse propagating 10 
nonlinear lengths, with an input pulse energy equal to that of the fundamental 
soliton, and with temporal durations ranging from 1 fs to 100 ps. The TOD 
has a value of (a) ps^/mm, (b) 0.1 ps^/mm, and (c) 0.2 ps^/mm. 

The analysis was then repeated for different ratio's of nonlinear lengths, and at numerous 
different TOD, ranging from to 0.2 ps^/mm. Again, the standard selected in order to 
consider a soliton pulse sustained within the waveguide was an output TBP deviation of 
less than 10% of the input TBP. As expected, with greater TOD, the minimum time has 
increased before the soliton has been disturbed. It was observed that the minimum time 
would increase almost linearly from no minimum time (as is the case of no TOD) up to 
100 fs, while the TOD increased linearly from to 0.15 ps^/mm. After this, increasing 
the minimum time appeared to plateau at 100 fs. While dramatic increases in the TOD 
have caused minimum times up to 500 fs, such dispersion was not observed in the practical 
PhCWG of interest. 

In the next step, a similar analysis was conducted relating to the disturbance of the 
phase of the output of the PhCWG. Throughout the study, a transform-limited pulse with 
no chirp was consistently used, where a phase of radian was applied to the entire pulse. 
On the output, it is expected that the phase would be shifted due to the SPM, and so the 
standard-deviation of the output pulse was determined. In the absence of TOD, there was 
no change in the standard-deviation of the output pulse, regardless of the input temporal 
duration; this is expected as the pulse has a consistent shape as a fundamental soliton. 
The standard selected for disturbance in the phase was a deviation of 0.1 radians of the 
standard-deviation of the output phase, when compared to the output phase at the longest 
100 ps temporal duration. By increasing the TOD, it was observed that the minimum time 
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Fig. 4. Data plots of soliton disturbance as a function of TOD for (a) 1 linear 
length, (b) 2 linear lengths, and (c) 10 linear lengths. Data plots of fundamental 
soliton phase disturbance as a function of TOD for (d) 1 linear length, (e) 2 
linear lengths, and (f) 10 linear lengths. 



would increase linearly with the TOD; the same as the minimum time for the soliton pulse 
stability. Unlike the pulse shape, however, the minimum time for phase stability would 
increase with increasing the nonlinear lengths. For a single nonlinear length, the minimum 
time was a maximum of 70 fs (for a TOD of 0.2 ps^/mm); whereas for ten nonlinear lengths, 
the minimum time for a stable phase was 1.3 ps (for the same TOD). It is observed from 
the numerical study that the soliton phase stability is much more sensitive to TOD than 
the soliton intensity shape. In practice, however, most PhCWG are less than a nonlinear 
length, so it is expected that there will be little phase disturbance as a result of TOD in 
practical applications of soliton pulse propagation within PhCWG. 



6. Conclusion 

In conclusion, a numerical study was conducted to theoretically understand what will 
happen to an ideal ultrafast fundamental soliton pulse propagating in a waveguide subjected 
to higher-order GVD. Simulations were conducted for fundamental solitons subjected to a 
large range of 3^^, 4^^, 5^^, and 6^^-order GVD, propagating with fs durations and 2^^-order 
GVD ranging from low-dispersion fibers to highly dispersive PhCWG. The study then 
turned to focus on practical applications of PhCWG, and it was determined that, when 
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the waveguide has TOD typical to a modern PhCWG, the sohton will become disturbed 
at pulses shorter than 100 fs. The phase of a fundamental soliton, which is subjected 
to self- focusing from SPM, was observed to become disturbed by the TOD at temporal 
durations ranging from 70 fs for a single nonlinear length, up to over a picosecond for 
ten nonlinear lengths; since most PhCWG are less than a full nonlinear length, TOD 
phase disturbance should not become an issue in practical applications for pulses under 
100 fs. In conclusion, the study demonstrated that fundamental soliton pulses longer than 
a few hundred fs could propagate undistorted in a practical PhCWG while subjected to 
higher-order GVD. By use of these PhCWG, one can accurately control and sustain these 
ultrafast, highly intense pulses over long distances, a phenomena that has a host of ap- 
plications ranging from low-energy optical data transfer to high-energy pulsed lasers systems. 
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